The measurement of properties of flowing vibrated yield stress fluids

ABSTRACT

A model of the viscoplastic boundary layer of a yield stress fluid is described and, based on which, there is provided a method comprises the steps of: vibrating a vibratory transducer at resonance at a first resonant mode in a yield stress fluid and making a first measurement of the resonant frequency; providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrating a vibratory transducer at resonance at the first resonant mode in the yield stress fluid and making a second measurement of the resonant frequency; and estimating the yield stress of the yield stress fluid based on the first and second measurements. The vibration to liquefy yield stress fluid around the one or more transducers may be provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement.

BACKGROUND

The present disclosure relates to techniques and apparatus for the measurement of fluid properties, and is concerned particularly with determining a yield stress of a non-Newtonian fluid.

In a Newtonian fluid, the shear stress is directly proportional to the shear rate, the constant of proportionality being the viscosity of the fluid. Therefore the viscosity as a single parameter can be used to model or define the relationship between shear stress and shear rate of the fluid, and therefore the flow behaviour of Newtonian fluids. Water is an example of a Newtonian fluid.

In a non-Newtonian fluid, the relationship between shear stress and shear rate is not so simple. The apparent viscosity of the fluid is found to vary with, for example, the shear stress or the shear rate. Fluids that exhibit non-Newtonian behaviour include tomato ketchup, mayonnaise and paint. The study of the flow of these types of fluids is the field of ‘rheology’.

Rheological models have been generated for non-Newtonian fluids, for which a small number of parameters can be obtained for a fluid, which may be used to determine the relationship between shear stress and shear rate in the fluid over a range of shear stresses or shear rates, or in other words the ‘apparent’ viscosity provided by the fluid under those particular conditions.

In the prior art, such parameters of rheological models are typically determined using a device that exerts shear forces in a fluid sample using a rotating or oscillating element. Complex viscoelastic parameters are usually determined by analysing the relationship between the shear rate and the developed shear stress. However, the devices use sensitive moving parts and their performance, and therefore the accuracy of measurement, can become affected by the environmental conditions in which they are used. Furthermore, the need to manage the fluid sample within a carefully defined volume, in order to calculate the shear rate accurately, generally makes these types of apparatus suitable only for the laboratory environment.

Of special interest is the real time measurement of the yield stress of a fluid. Yield stress is defined as the stress required in order for solid matter to flow and represents the point of plastic deformation of the solid. It is a feature of many fluids and has particular importance in the manufacture of foods, paints and petro-geological fluids, for example, but also in naturally occurring materials such as muds and sediments.

The details and origins of yield stress behaviour of fluids is a subject of interest to researchers. Traditional approaches to determining a yield stress are set out in “Yield Stress in Foods: Measurements and Applications” of Sun et al (International Journal of Food Properties, volume 12, pp 70-101, 2009), and in “Understanding Yield Stress Measurements”, a white paper produced by Malvern Instruments Limited of Grovewood Road, Malvern, Worcestershire, UK.

Traditional approaches measure yield stress of a static sample.

For example, in the extrapolation approach, there is extrapolation of the shear stress versus shear rate data obtained from conventional rheometers. The experimental data is extrapolated back to zero shear rate to obtain the yield stress value at the shear stress intercept.

For example, in the stress relaxation approach, the fluid material is first sheared at either constant shear rate or constant shear stress in a conventional rotational viscometer, followed by bringing the material to rest either gradually or suddenly. The yield stress is then measured as residual stress remaining in the fluid upon cessation of flow.

For example, in the creep/recovery response approach, a constant shear stress is applied in steps to the fluid material. If stresses applied are below the yield stress, the test material behaves as an elastic solid with a complete recovery upon removal of stress and it will not flow. The yielding point can be detected from a drastic change of slope in the time vs. shear strain curve.

For example, in the shear stress ramp approach, a gradual step increase is applied to the sample. The instantaneous (or apparent) viscosity is monitored for the presence of an inflexion point, which indicates onset of flow, and the yield stress.

For example, in the cone penetrometer approach, a metal cone of specific dimensions is forced into the surface of test specimen. The cone gradually slow down until it comes to rest. The yield stress is calculated as a function of the penetration depth when the cone stops.

For example, in the dynamic oscillation approach, the material is subjected to a sinusoidal strain and the resulting stress is measured as a function of both time and frequency. Depending on the strain amplitude, small deformations occur within the material. The presence of a low frequency plateau in the curve can be correlated to the yield stress.

For example, in the inclined plane approach, a uniform fluid layer is placed on an initially horizontal plane, followed by progressively increasing the angle of inclination of the plane until a critical value is reached, whereby the fluid starts to flow. The yield stress is correlated with the angle of inclination for the fluid to flow.

For example, in the vane technique, a vaned device is immersed in the material. The vane geometry consists of a vane spindle with typically four to eight thin blades arranged at equal angles, cantered on a narrow cylindrical shaft of a rheometer or viscometer. The vane can be operated under controlled shear-rate mode or controlled shear stress mode; in both cases the yield stress is the minimum stress required for continuous rotation of the vane.

For example, in the squeezing flow approach, a sample of fluid is placed between the two parallel test surfaces of a bench rheometer or viscometer and squeezed in either a ‘constant volume’ or ‘constant area’ arrangement. The determination of the yield stress from squeezing flow is based on the residual stress after shear and relaxation.

For example, in the plate method, a plate is immersed into a container filled with sample material. An attached balance records the force required to create plate movement and this is combined with the buoyancy force to determine yield stress.

For example, in the slump test, a cylindrical mould is filled with the test fluid and the mould is lifted off to allow the material to collapse and flow under its own weight. The slump height, which is the difference between the initial and final heights is measured and provides an inverse relationship with yield stress.

Unlike the above-described techniques, another traditional approach allows for yield stress measurement of a flowing sample, by way of the use of magnetic resonance imaging (MRI), whereby a material sample flowing at a constant velocity through a tube is flow imaged using an MRI scanner. The MRI scanner is able to measure the flow profile (flow velocity) of the fluid across the pipe. The flow profile is combined with an estimate of stress from the differential pressure along the pipe to provide a measurement of yield stress.

The disadvantage of almost all of the above traditional approaches to determining the yield stress of a fluid is the requirement for a captive sample in a static (non-flowing) environment. This makes these methods unsuitable for real-time measurement. The traditional methods requiring captive samples are largely confined to intermittent spot measurements and the laboratory environment.

It would be desirable to provide fast real-time rheology measurement of yield stress fluids under live conditions. This would allow detailed assessment of large material volumes and improved capability for their control and management. In particular, in process industries, in-line measurement is advantageous in that it allows for continuous product control. A feedback loop may allow parameters of the process to be adjusted based on in-line measurements of yield stress. While the traditional approaches using captive static samples may be used to take periodic measurements, there is inevitably a measurement lag. The provision of faster measurements in-line would reduce or eliminate the measurement lag and thus improving the controllability of the process, which can increase product consistency and reduce waste.

In other applications such as the measurement of rheological properties of drilling muds and lubricants, the fast real time rheological measurement of yield stress under live conditions may improve efficiency by allowing faster intervention to correct drift in fluid properties away from desired or permitted ranges.

In marine, civil engineering, mining and geotechnical fields it would be desirable to provide fast measurement of flowing materials with important rheology features such as poured concretes and cements, slurries and muds, in situ, without taking captive static samples, providing increased accuracy due to ease of taking multiple data points, while reducing or eliminating delays involved with spot testing and laboratory turnaround times.

SUMMARY

According to a first aspect of the invention, there is provided a method of estimating the yield stress of a yield stress fluid using one or more vibratory transducers, the method comprising: vibrating a vibratory transducer in the yield stress fluid and making a first measurement; providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrating a vibratory transducer in the yield stress fluid and making a second measurement; estimating the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency. Vibration sources of an appropriate amplitude and frequency may be applied locally to the yield stress material surrounding the viscometer transducer to create liquefaction. This may be performed in flowing or static conditions. Yield stress is estimated by either the reduction in the apparent viscosity or a change in frequency due to a change in mass loading at the transducer surface. The vibration source may be the resonance of the transducer itself as used to make a viscosity measurement and in the context of the present disclosure, this is termed ‘intrinsic stimulation’. Alternatively it may be a separate vibration. The combination of high frequency and high amplitude of vibration may cause local liquefaction at and around the transducer surface, detected as a viscosity and frequency change.

Preferably, the vibration to liquefy a yield stress fluid around the one or more transducers is provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement.

More preferably, the yield stress is estimated based on the difference between the first and second measurements, the difference scaled by the difference in amplitude between the first and second measurements.

Alternatively, method further comprises obtaining third and fourth measurements at a lower frequency and without liquefaction, wherein the yield stress is estimated based on the difference between the first measurement scaled by the third measurement and the second measurement scaled by the fourth measurement, said difference being scaled by the difference in amplitude between the first and second measurements. The use of measurements obtained at the lower frequency without liquefaction allows temperature effects and other common mode fluctuations to be ratiometrically reduced or eliminated. The first, second, third and fourth measurements may be made using a single multi-mode vibratory transducer or using multiple transducers. For example, the first and second measurements may be made using a first vibratory transducer and the third and fourth measurements may be made using a second vibratory transducer.

In a preferred embodiment, a method of estimating the yield stress of a yield stress fluid using one or more vibratory transducers comprises the steps of: vibrating a vibratory transducer at resonance at a first resonant mode in the yield stress fluid and making a first measurement of the resonant frequency; providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrating a vibratory transducer at resonance at the first resonant mode in the yield stress fluid and making a second measurement of the resonant frequency; estimating the yield stress of the yield stress fluid based on the first and second measurements. Preferably, the vibration to liquefy yield stress fluid around the one or more transducers is provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement. More preferably, the yield stress is estimated based on the difference between the first and second measurements of the resonant frequency scaled by the difference in amplitude between the first and second measurements. Alternatively, the method may further comprise: vibrating a vibratory transducer at resonance at a second resonant mode in the yield stress fluid and making a third measurement of the resonant frequency, the second resonant mode corresponding to a lower resonant frequency than the first resonant mode; and vibrating a vibratory transducer at resonance at the second resonant mode in the yield stress fluid and making a fourth measurement of the resonant frequency, wherein the yield stress is estimated based on difference between the first frequency measurement scaled by the third frequency measurement and the second frequency measurement scaled by the fourth frequency measurement, scaled by the difference in amplitude between the first and second measurements. Advantageously, the use of the measurements of the second resonant mode allows temperature effects and other common mode fluctuations to be ratiometrically reduced or eliminated.

In another preferred embodiment, a method of estimating the yield stress of a yield stress fluid using one or more vibratory transducers comprises the steps of vibrating a vibratory transducer in the yield stress fluid at a first frequency and making a first measurement of the degree of damping; providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrating a vibratory transducer the first frequency and making a second measurement of the degree of damping; estimating the yield stress of the yield stress fluid based on the first and second measurements of degree of damping. Preferably, the vibration to liquefy yield stress fluid around the one or more transducers is provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement. More preferably, the yield stress is estimated based on the difference between the first and second measurements of the degree of damping scaled by the difference in amplitude between the first and second measurements. Alternatively, the method may further comprise: vibrating a vibratory transducer at a second frequency in the yield stress fluid and making a third measurement of the degree of damping, the second frequency being lower than the first frequency; and vibrating a vibratory transducer at the second frequency in the yield stress fluid and making a fourth measurement of the degree of damping, wherein the yield stress is estimated based on difference between the first degree of damping measurement scaled by the third degree of damping measurement and the second viscosity measurement scaled by the fourth degree of damping measurement, scaled by the difference in amplitude between the first and second measurements.

Preferably, the measurements of the degree of damping of the yield stress fluid comprise measurements of the viscosity of the yield stress fluid.

According to a further aspect of the invention, there is provided an apparatus for estimating the yield stress of a yield stress fluid using one or more vibratory transducers, the apparatus comprising: one or more vibratory transducers; a processing module configured to: vibrate a vibratory transducer at resonance at a first resonant mode in the yield stress fluid and make a first measurement of the resonant frequency; provide a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrate a vibratory transducer at resonance at the first resonant mode in the yield stress fluid and make a second measurement of the resonant frequency; estimate the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency.

According to a further aspect of the invention, there is provided an apparatus for estimating the yield stress of a yield stress fluid using one or more vibratory transducers, the apparatus comprising: means for vibrating a vibratory transducer in the yield stress fluid and making a first measurement; means for providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; means for while said portion of the yield stress material is liquefied, vibrating a vibratory transducer in the yield stress fluid and making a second measurement; means for estimating the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency.

According to a further aspect of the invention, there is provided a non-transitory computer-readable medium having stored thereon instructions that, when executed by one or more processors, cause the one or more processors to: vibrate a vibratory transducer in the yield stress fluid and making a first measurement; provide a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrate a vibratory transducer in the yield stress fluid and make a second measurement; estimate the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency.

In preferred embodiments, the apparatus comprises means for carrying out one or more of the above-described methods or indeed any of the methods described herein.

According to a further aspect, there is provided a non-transitory computer-readable medium having stored thereon instructions that, when executed by one or more processors, cause the one or more processors to: vibrate a vibratory transducer in the yield stress fluid at a first frequency and make first measurement of the viscosity; provide a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrate a vibratory transducer the first frequency and make a second measurement of the viscosity; estimate the yield stress of the yield stress fluid based on the first and second measurements of viscosity.

Aspects of the invention relate to the measurement of one or more of the following properties of Newtonian and non-Newtonian fluids: yield stress, viscosity (at one or more equivalent shear rates), power law index n, consistency factor K, flow velocity, and the detection of onset of flow from initially static conditions. In addition, aspects of the invention may be applicable in both static and flowing conditions.

Aspects of the invention may be used in a wide range of settings, for example: as an in-line harsh process-tolerant device in pipes, vessels, open channels; for the taking of field measurements in marine and earth science environments; in the laboratory with continuous, automated or discrete sample measurements; as a portable device for both laboratory, field and at-line use; and as a component of a system incorporating measurement of these fluid properties including resonators based on micro-electro-mechanical systems (MEMS) and nano-electro-mechanical systems (NEMS) scale architectures.

Aspects of the invention are based on the use of one or more viscosity measurement transducers and the modulation of the rheological state of the fluid at the transducer measurement surface through one or more of the formation of a boundary layer as a result of fluid flow and the liquefaction of non-yielded solid material by vibrational agitation.

In some cases, aspects of the invention can be performed using a conventional viscometer design, such as rotational viscometers including rotational viscometers of the cone and plate type and the rotational cylinder type. Other conventional viscometer designs with which aspects of the invention may be performed include differential pressure viscometers and falling-ball-type viscometers. However, compared with traditional approaches to measuring yield stress, aspects of the invention may have particular application for the surface-loaded attributes of resonant viscometer transducers, especially such types which work by the formation and dissipation of a pure shear wave.

Compared with traditional approaches to measuring yield stress that for the most part require a captive static sample of fluid, aspects of the invention may have particular application in the continuous rheological measurement of materials in the flowing condition, while also being capable of use with static materials.

In particular, aspects of the invention may have application in the monitoring of yield stress or other fluid properties in-line in a process in which the yield stress and its variation is a value of interest, and further substances may be added in greater or lesser amounts according to the estimated yield stress of the fluid and any particular desired yield stress or other fluid property of the yield stress fluid.

BRIEF DESCRIPTION OF THE DRAWINGS

Aspects of the invention will be described in more detail by way of example only with reference to the accompanying drawings. Components within the drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating principles.

FIG. 1 is a cross-sectional view of a yield stress fluid in a pipe, wherein the fluid is at rest;

FIG. 2 is a cross-sectional view of a yield stress fluid in a pipe, wherein the fluid is flowing;

FIG. 3 is a cross-sectional view of a static yield stress fluid in which there is located a viscosity transducer in the form of a bob attached to a shaft, the transducer being operable at at least two frequencies, f1, f2;

FIG. 4 is a cross-sectional view of a flowing yield stress fluid in which there is located a viscosity transducer in the form of a bob attached to a shaft, the transducer being operable at at least two frequencies, f1, f2;

FIG. 5 is a cross-sectional view of a flowing yield stress fluid in which there is located a first viscosity transducer in the form of a bob attached to a shaft, the transducer being operable at a first frequencies, f1, and a second viscosity transducer in the form of a bob attached to a shaft, the transducer being operable at a second frequencies, f2;

FIG. 6 is a cross-sectional view of a flowing yield stress fluid in which there is located a multiple-frequency rod resonator capable of vibrating in torsional, lateral and longitudinal modes;

FIG. 7 is a cross-sectional view of a flowing yield stress fluid in which there is located a multiple-frequency disc resonator capable of vibrating in torsional, lateral and longitudinal modes, the axis of the disc being perpendicular to the flow direction;

FIG. 8 is a cross-sectional view of a flowing yield stress fluid in a vibrating tube viscometer that is capable of vibrating in torsional, lateral and longitudinal modes;

FIG. 9 is a cross-sectional view of a yield stress fluid a vibrating vessel viscometer that is capable of vibrating in torsional, lateral and longitudinal modes;

FIG. 10 is a cross-sectional view of a static yield stress fluid in which there is located a viscosity transducer in the form of a bob attached to a shaft, wherein no boundary later forms;

FIG. 11 is a cross-sectional view of a flowing yield stress fluid in which there is located a viscosity transducer in the form of a bob attached to a shaft, wherein a boundary later forms around the viscosity transducer;

FIG. 12 illustrates the tri-viscous model by showing a yield stress fluid flowing past a multi-frequency resonator surface and forming a boundary layer, in which shear forces due to the flow of yield stress material break up the solid plug, forming a graded zone from solid to liquid across the boundary layer;

FIG. 13 shows, a wave emanating from the surface of the transducer that will firstly traverse the liquid region within the boundary layer, Ld, and therefore experience the liquid viscosity, V_(L). Where the wave propagation depth is sufficiently long, the wave may reach the transitional region and thereby register the higher viscosity V_(T) of the material in the transitional region;

FIG. 13, a wave emanating from the surface of the transducer will firstly traverse the liquid region within the boundary layer, Ld, and therefore experience the liquid viscosity, V_(L). Where the wave propagation depth is sufficiently long, the wave may reach the transitional region and thereby register the higher viscosity V_(T) of the material in the transitional region;

FIG. 14 shows an alternative configuration in which a third wave is generated by way of a third measurement at a third frequency;

FIG. 15 shows flow modulation by varying the viscometer body orientation to the flow;

FIG. 16 shows flow modulation by viscometer position relative to a nearby surface. Either or both of the surface and transducer can move;

FIG. 17 shows flow modulation by a local flow obstructer;

FIG. 18 shows flow modulation by modulating a container position;

FIG. 19 illustrates liquefaction by intrinsic vibration, with example lateral, longitudinal and torsional vibrational modes shown;

FIG. 20 illustrates liquefaction by extrinsic vibration, the vibration being provided at a vessel or conduit wall, with example lateral, longitudinal and torsional vibrational modes shown;

FIG. 21 illustrates liquefaction by extrinsic vibration, whereby vibration to cause liquefaction is provided by vibration sources located elsewhere in the yield stress material;

FIG. 22 shows schematically an apparatus for carrying out one or more techniques of the present disclosure;

FIG. 23 shows schematically an apparatus for carrying out one or more techniques of the present disclosure;

FIG. 24 shows schematically an apparatus for carrying out one or more techniques of the present disclosure;

FIG. 25 shows schematically an apparatus for carrying out one or more techniques of the present disclosure;

FIG. 26 shows schematically an apparatus for carrying out one or more techniques of the present disclosure;

FIG. 27 shows schematically an apparatus for carrying out one or more techniques of the present disclosure, wherein a viscosity transducer is pivotable from a first position perpendicular to the flow direction, to a second position rotated toward the flow direction;

FIG. 28 shows a transducer with a smooth profile in open flow of yield stress material;

FIG. 29 shows the transducer of FIG. 28 in a flow of yield stress material through a pipe or vessel;

FIG. 30 shows a transducer with a contoured profile in an open flow of yield stress material;

FIG. 31 shows the transducer of FIG. 30 with a contoured profile in a flow of yield stress material through a pipe or vessel;

FIG. 32 shows the smooth transducer of FIG. 28 in a flow of yield stress material through a pipe or vessel, the walls of the pipe or vessel being axially aligned with the flow direction and the shaft of the transducer, wherein recesses are provided in the walls of the pipe or vessel;

FIG. 33 shows the transducer of FIG. 30 with a contoured profile in a flow of yield stress material through a pipe or vessel, wherein a further boundary layer develops at the walls of the pipe or vessel, the walls of the pipe or vessel being axially aligned with the flow direction and the shaft of the transducer, wherein recesses are provided in the walls of the pipe or vessel;

FIG. 34 shows the transducer of FIG. 28 in open flow of yield stress material but, unlike in FIG. 28, the transducer is aligned perpendicular to the flow direction;

FIG. 35 shows the transducer of FIG. 34 aligned perpendicular to flow in a pipe, wherein the transducer is partially retracted into a recess of the pipe;

FIG. 36 shows a resonant disc transducer partially retracted into a recess;

FIG. 37 shows a series of side views of five bob profiles as possible designs for a transducer of the shaft-and-bob type;

FIG. 38 shows a side view of two further bob designs for a transducer;

FIG. 39 shows two further bob designs for a transducer in perspective and cross-sectional view;

FIG. 40 shows a further bob design for a transducer in perspective and cross-sectional view.

DETAILED DESCRIPTION

A yield stress fluid effectively behaves like a solid when shear stresses in the fluid are below the yield stress. This happens, for example, when the fluid is at rest or when a volume of fluid is moving with uniform velocity. When a flowing material possessing yield stress encounters a surface, a velocity gradient is developed in the region of material near the surface. This gradient leads to the formation of shear stresses within the fluid. Liquefied regions will form where these stresses exceed the yield stress of the material. The liquefied layer regions developing as a result of the velocity gradient in the material near the surface is the boundary layer (BL).

FIG. 1 shows a yield stress fluid in a pipe. There is no flow and so the fluid is at rest. Therefore there is no velocity gradient in the fluid and so it is unyielded, i.e. solid, at all locations.

FIG. 2 shows a yield stress fluid undergoing flow in a pipe, where the solid material shears at the pipe wall forming a liquid layer. Due to friction and the velocity of the fluid relative to the pipe wall, there is a velocity gradient leading to yielded material, i.e. liquid, in a boundary layer adjacent the pipe wall. Further from the pipe wall, the velocity gradient decreases leading to a reduction in shear stress to below the yield stress, causing the flowing material to solidify. This flow is known as ‘plug flow’.

The unyielded material is considered to be in the solid regime. The yielded material is considered to be in the liquid regime and to have a specific depth.

There is a reduction of apparent viscosity as the fluid liquefies. The transition from liquid to solid at the periphery of the boundary layer is complex and the subject of ongoing research. Material that may assist in understanding the present disclosure and the complex behaviour in the boundary layer of a yield stress material undergoing flow includes Boujlel et al “Boundary Layer In Pastes—Displacement Of A Long Object Through A Yield Stress Fluid” (Journal of Rheology, volume 56, 2012, doi:10.1122/1.4720387).

Without wishing to be bound by theory, it is believed that part or even all of the liquid layer could be considered a region of graded viscosity increasing rapidly into the solid zone.

The depth of the boundary layer, d, is found to be a function of flow velocity U but may also be affected by viscosity and yield stress.

Immersing wholly or at least partially a vibratory viscosity transducer into a fluid can enable viscosity and resonant frequency to be measured in real time.

FIGS. 3, 4 and 5 show by way of example a viscosity transducer in the form of a bob attached to a shaft, the arrangement located in a yield stress material and vibrating torsionally relative to a central longitudinal axis.

FIG. 3 shows that, with the material at rest, the environment around the transducer is such that there is no boundary layer—the material is wholly solid.

FIG. 4 shows that material undergoing flow will shear against the transducer surface causing the formation of a liquefied boundary layer.

The transducer is capable of operating at least two frequencies. A first frequency of operation may be relatively low, such as below 400 Hz. A second frequency of operation may be relatively high, such as above 1500 Hz. It is to be noted that these frequencies are by way of example only and that the invention may be put into effect using frequency ranges other than the examples set out here.

Vibration at each frequency produces a shear wave that will propagate in the boundary layer.

The higher frequency wave has a relatively short penetration depth, typically covering the most highly sheared liquefied region and relatively little, if any, of the solid regime. This results in a lower viscosity measurement and a lower mass loading at the transducer.

The lower frequency wave has a relatively long penetration depth, typically covering the most highly sheared liquefied region but also a greater portion of the solid regime, the particular amount being influenced by the boundary layer depth. Depending on the yield stress, the fluid will appear more solid at the lower frequency, resulting in a higher viscosity measurement and increased mass loading at the transducer.

By cancelling out the effect of the liquefied region to leave only the effect of the solid region that is influenced by boundary layer depth, an estimate of the yield stress can be made. In particular, an estimate of the yield stress can be made as a differential function of the viscous loss readings at the two frequencies.

Advantageously, this approach does not require many of the assumptions made by other approaches to determining a yield stress. For example, it is not imperative that the fluid obeys the Cox-Merz rule, whereby a steady-state shear viscosity at a given shear rate is approximated by the dynamic viscosity at the same frequency. The Cox-Merz rule is an empirical rule and applies only for certain non-Newtonian fluids. Due to not being reliant on the Cox-Merz rule, this approach has a wider application and is not limited to fluids for which the Cox-Merz rule is a valid approximation.

In addition, compared with existing approaches in which parameters are fit to fluid models such as the Herschel-Bulkley or Casson models, the yield stress representing the zero-shear-rate intercept, this approach reduces or eliminates variation due to thermal effects, noise, unsteady flow rate. This is consequence of the form of the solution as a differential function of the viscous loss readings—by subtracting one measurement from another, common-mode errors or variations are reduced or eliminated.

Viscosity of a yield stress fluid may be measured at two different vibrational (resonant) frequencies. For example, the two frequencies may be a relatively low frequency of 400 Hz and a relatively high frequency of 1500 Hz. This may be achieved with two single frequency viscometers in relatively close proximity to each other as shown in FIG. 5 or alternatively as a single device operating as a dual-frequency viscometer as shown in FIG. 4.

The viscosity transducer is not limited to the form shown in FIG. 4. In this example, the viscosity transducer comprises a bob attached to a shaft. The bob may take many geometric forms including, but not limited to, a cylinder, disc, or sphere.

FIG. 6 shows a multiple-frequency rod resonator capable of vibrating in torsional, lateral and longitudinal modes. Such a device may act as a dual-frequency viscometer by selecting two modes having different resonant frequencies.

FIG. 7 shows a multiple-frequency disc resonator capable of vibrating in torsional, lateral and longitudinal modes. Such a device may act as a dual-frequency viscometer by selecting two modes having different resonant frequencies.

In such cases the viscosity transducers are preferably configured to vibrate in a torsional mode but lateral and longitudinal modes of vibration can also be used.

Alternatively the transducer may be the actual fluid container as in the case of a resonant tube device or a vibrating vessel of any shape.

FIG. 8 shows a vibrating tube viscometer capable of vibrating in torsional, lateral and longitudinal modes.

FIG. 9 shows a vibrating vessel viscometer capable of vibrating in torsional, lateral and longitudinal modes.

As before, the viscosity transducers are preferably configured to vibrate in a torsional mode but lateral and longitudinal modes of vibration can also be used.

Using the bob-and-shaft-type resonator as an example, but applicable to all cases: in the absence of an additional liquefaction vibration source, under non-flowing conditions the transducer will effectively detect a solid material as shown in FIG. 10 and, when the material flows, a liquid boundary layer region will form as shown in FIG. 11.

In the drawings of FIGS. 3 to 11, the boundary layer is seen as a clearly marked zone (broadly spaced hatching aligned top-left to bottom-right) against the solid material (more narrowly spaced hatching aligned top-right to bottom-left).

However, without wishing to be bound by theory, it is believed that the boundary layer is more complex than this, and particularly that it is a region where a highly liquid region graduates towards a solid region, with a mixture of both liquid and solid elements in a transitional region lying between the solid and highly liquid regions.

Some previously published models refer to a two-state nature of the boundary layer as ‘bi-viscous’. Once again, without wishing to be bound by theory, it is believed that the change from liquid to solid is better modelled by three separate viscosity zones, liquid, transitional and solid, because it is believed that bi-viscous models do not easily accommodate a non-sudden change of state from liquid to solid across the boundary layer.

The use of three separate viscosity zones is presented for the first time in this disclosure and is termed a ‘tri-viscous’ model of the boundary layer.

FIG. 12 illustrates the tri-viscous model by showing a yield stress fluid flowing past a multi-frequency resonator surface and forming a boundary layer, in which shear forces due to the flow of yield stress material break up the solid plug, forming a graded zone from solid to liquid across the boundary layer. Three shear waves are shown. The lowest frequency wave of angular frequency ω1 penetrates furthest into the boundary layer. The highest frequency wave of angular frequency ω3 remains close to the resonator surface. The intermediate frequency wave of angular frequency ω2 penetrates to a distance somewhere in-between the lowest and highest frequency waves.

With reference to FIG. 12, the bulk of the material flows with a velocity U₀ and a velocity profile is formed, the velocity local to the surface having a maximum value of U₀ outside the boundary layer and falling to zero at the resonator surface.

The velocity profile gives rise to a shear rate profile, with the shear rate varying from a maximum at the surface to zero outside the boundary layer.

At very low or zero shear rate the material remains a solid with an apparent viscosity, V_(S).

As the shear rate increases, some of the solid material yields, leading to a fragmented dispersion of solid material. This fragmented dispersion of solid material frustrates the propagation of the shear wave leading to a high apparent fluid viscosity, V_(F), as shown in the viscosity profile of viscosity vs distance from the transducer surface. This leads to a transitional viscosity region, V_(T).

At higher shear rates, more of the solid material yields causing a reduction in the fragmented dispersion of solid material and the material increasingly behaves as a liquid, easing wave propagation and leading to a reduction in apparent viscosity V_(F) to the apparent liquid viscosity, V_(L).

The increased apparent viscosity of V_(T) is a result of the corruption of the homogeneity of the liquid layer by fragmentation nearing the boundary layer. V_(T) has a higher value than V_(L), or even V_(S), as the propagation path of the wave through the fragments of solid material that are dispersed in liquid is more dissipative than through a homogeneous liquid or solid.

FIG. 13 shows the tri-viscous model more clearly, without illustrating the fragmentary nature of the transitional region.

The tri-viscous model can be used to estimate the yield stress of a material.

According the Herschel-Bulkley model,

σ=σ₀ +Kγ ^(n),  (1)

in which σ represents the shear stress, σ₀ represents the yield stress, K is a fluid-dependent parameter termed the ‘consistency’, γ represents the shear rate, and n is a fluid-dependent parameter termed the ‘power law index’.

Dividing equation (1) by the shear rate γ results in:

$\begin{matrix} {{\frac{\sigma}{\gamma} = {\frac{\sigma_{0}}{\gamma} + {K\; \gamma^{n - 1}}}},} & (2) \end{matrix}$

Looking at the sources of the shear rate, it has two components:

γ=γ_(F)+γ_(ω),  (3)

where γ_(F) is the shear rate due to the bulk fluid flow and γ_(ω) is the shear rate due to vibration of the resonator.

By definition, the apparent fluid viscosity V_(F) is the ratio of shear stress and shear rate. Writing equation (2) in terms of V_(F): results in

V _(F) =V _(YS) +V _(L),  (4)

where a portion V_(YS) of the apparent viscosity V_(F) is due to the yield stress and a portion V_(L) is due to the liquefied material. In the case of viscosity due to yield stress, the viscous wave dissipation is a function of the inhomogeneity caused by shear flow, and so γ=γ_(F)

$\begin{matrix} {{V_{YS} = {\frac{\sigma_{0}}{\gamma} = \frac{\sigma_{0}}{\gamma_{F}}}},} & (5) \end{matrix}$

In the case of resonant viscosity measurements, γ_(ω) is significantly greater than γ_(F) and so γ→γ_(F). Applying the Cox-Merz rule, γ→γ_(F)→ω. Therefore the apparent viscosity of liquefied material V_(L) is given by

V _(L) =Kω ^(n-1)  (6)

Combining equations 4, 5, and 6 results in an expression for the apparent fluid viscosity in terms of the yield stress, shear rate due to flow, consistency, frequency and power law index:

$\begin{matrix} {V_{F} = {\frac{\sigma_{0}}{\gamma_{F}} + {K\; \omega^{n - 1}}}} & (7) \end{matrix}$

The shear rate due to flow γ_(F). is a function of depth and other parameters, and can be defined in terms of a shear rate at the surface, γ₀, and a depth function ƒ(d/B_(d)) as

γ_(F)=γ₀/ƒ(d/B _(d))  (8)

where the depth function is a somewhat arbitrary function developed to provide an expected variation in shear rate across the boundary layer, given by:

ƒ(d/B _(d))=(1−(d/B _(d))^(R))⁻¹  (9)

where d is the wave propagation depth, Bd is the boundary layer depth, and R is a fluid-dependent shear rate index. This leads to

$\begin{matrix} {V_{F} = {{\frac{\sigma_{0}}{\gamma_{0}}{f\left( {d\text{/}B_{d}} \right)}} + {{V_{L}(\omega)}.}}} & (10) \end{matrix}$

A shear wave propagating from the resonant viscometer will dissipate to a skin depth a distance d from the surface. The propagation depth in a viscoelastic fluid varies with the loss tangent tan(δ) of the material. But, for illustration, the simpler Newtonian interpretation of this distance may be used, in which the skin depth is given by:

$\begin{matrix} {d = \sqrt{\frac{2V}{\omega\rho}}} & (11) \end{matrix}$

where ρ is the fluid density.

According to equation (11), within a given viscosity-density regime, the propagation depth can be assumed relatively constant for a given frequency.

With reference to FIG. 13, a wave emanating from the surface of the transducer will firstly traverse the liquid region within the boundary layer, Ld, and therefore experience the liquid viscosity, V_(L). Where the wave propagation depth is sufficiently long, the wave may reach the transitional region and thereby register the higher viscosity V_(T) of the material in the transitional region.

Two waves may be configured to emanate from the surface of the transducer at different frequencies, ω1 and ω2. The relatively low frequency wave at angular frequency ω1 propagates further than the relatively high frequency wave at angular frequency ω2. The frequencies are selected such that the low frequency wave experiences the viscosity of the transitional region V_(T) at a penetration distance d₁ that is greater than the depth Ld of the liquid region in the boundary layer and the high frequency wave has a short penetration depth d₂ and typically covers only the highly sheared, liquefied region in the boundary layer.

From equation (10),

$\begin{matrix} {{V\; 1} = {V_{F} = {{\frac{\sigma_{0}}{\gamma_{0}}{f\left( {d_{1}\text{/}B_{d}} \right)}} + {{V_{L}({\omega 1})}.}}}} & (12) \end{matrix}$

In normal conditions, d (wave propagation depth), Bd (boundary layer depth), and R (fluid-dependent shear rate index) and γ₀ (shear rate at surface) can be considered relatively constant, and so equation (12) simplifies to

V1=ασ₀ +V _(L)(ω1)  (13)

for some empirical value α.

In the case of the high frequency wave, the measured viscosity V2 will mostly represent the liquid-only viscosity, according to which equation (10) simplifies to:

V2=V _(L)(ω2).  (14)

In the liquid region, the viscosity will vary with shear rate and therefore the liquid-region viscosity of the low frequency wave V_(L)(ω1) will not equal liquid-region viscosity of the high frequency wave V_(L)(ω2). But, V_(L)(ω1) can be related to V_(L)(ω2) via the power law equation using the power law index n, according to which:

$\begin{matrix} {{V_{L}({\omega 1})} = {{V_{L}({\omega 2})}{\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}.}}} & (15) \end{matrix}$

Substituting equation (14) into equation (15) leads to

$\begin{matrix} {{V_{L}({\omega 1})} = {V\; 2{\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}.}}} & (16) \end{matrix}$

Substituting equation (16) into equation (13) leads to

$\begin{matrix} {{V\; 1} = {{\alpha\sigma}_{0} + {V\; 2{\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}.}}}} & (17) \end{matrix}$

Rearranging leads to an expression for yield stress of the fluid given by

$\begin{matrix} {{\sigma_{0} = {\alpha^{\prime}\left( {{V\; 1} - {V\; 2\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}}} \right)}},} & (18) \end{matrix}$

where V1 is the measured viscosity at frequency ω1, V2 is the measured viscosity at frequency ω2, n is the power law index and a′ is an empirical scaling constant determined by, for example, calibration with a test sample.

For known and fixed frequencies ω1 and ω2, the yield stress can even be written as the linear combination of V1 and V2 according to the expression

σ₀ =α′V1+a″V2  (19)

in which α′ and α″ are each empirical values. In fact, equation (19) is a generalized case of equation (18).

The ratio of α″ to α′ is given by the expression:

$\begin{matrix} {- {\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}.}} & (20) \end{matrix}$

Estimating the yield stress of a material using equation (18) requires a value to be supplied for the power law index n of the fluid. For many materials, predetermined values of the power law index n are available. For example, n is usually taken to be 0.5 for molten chocolate.

FIG. 14 shows an alternative configuration in which the value of n can be obtained as a real-time measurement by using a viscosity measurement at a third frequency ω3 that is higher than ω2. If ω2 is 1500 Hz, ω3 might be 5 kHz, say.

Both ω2 and ω3 are chosen to have sufficiently high frequency values that they will have a propagation depth lower than the liquid depth Ld. The propagating waves will then only experience the liquefied state within the layer. Therefore V2=V_(L)(ω2) and V3=V_(L)(ω2).

For the purpose of determining n, the yield stress can be ignored and we use the power law expression in equitation (15) to obtain a relationship between V2, V3, ω2, ω3, and n:

$\begin{matrix} {{{V\; 2} = {V\; 3\left( \frac{\omega 2}{\omega 3} \right)^{n - 1}}},} & (21) \end{matrix}$

which can be rearranged to obtain n by the expression:

$\begin{matrix} {n = {\frac{\log \left( {V\; 2\text{/}V\; 3} \right)}{\log \left( {{\omega 2}\text{/}{\omega 3}} \right)} + 1}} & (22) \end{matrix}$

Once a value of n is obtained, it can be applied in, for example, equation (18) to estimate the yield stress.

The shear rate at the surface γ₀ has been found to vary as a function of flow velocity U₀ and boundary layer depth B_(d) as follows:

$\begin{matrix} {\gamma_{0} = \frac{\beta^{\prime} \cdot U_{0}}{B_{d}}} & (23) \end{matrix}$

where β′ is an empirical constant

Inserting equation (23) into equation (12), V1 can be seen as a function of flow velocity U₀ by way of the following expression:

V1=V _(F)=σ₀(B _(d) /β′U ₀)ƒ(d ₁ /B _(d))+V _(L)(ω1).  (23A)

According to this expression, V1 increases with decreasing flow rate. By contrast, V2, in equation (14) is not seen to vary with flow rate.

Therefore, by dividing one viscosity by the other. i.e. by obtaining, VR=V2/V1, an estimate of flow velocity, U₀ is obtained.

$\begin{matrix} {U_{0} = \frac{{VR} \cdot {kA}}{{kB} - {{kC} \cdot {VR}}}} & (24) \end{matrix}$

where kA, kB and kC are empirically determined calibration parameters.

While the expression for determining a yield stress in equation (18) provides advantages in that it is a closed-form analytical expression that allows the value to be computed directly, this is not the only possible approach to calculating yield stress.

For example, equation (1) representing the Herschel-Bulkley model can be rewritten as

σ₀ =σ−Kγ ^(n),  (25)

The shear stress σ for a given shear rate γ is a function of viscosity as follows:

σ=V _(L)·γ,  (26)

Shear rate can be derived from the viscometer frequency of vibration using the Cox-Merz rule, where shear rate is directly equivalent to the angular frequency ω. Thus combining equations (25) and (26) results in:

σ₀ =V _(L) ·ω−Kω ^(n)  (27)

Taking viscosity and frequency readings from two resonant viscometers or a single multi-mode device leads to the simultaneous equations (assuming the power law index n is known):

σ₀ =V _(L)(ω1)·ω1−Kω1^(n)

σ₀ =V _(L)(ω2)·ω2−Kω2^(n)  (28)

These equations can be solved to find the yield stress σ₀ and the consistency K. If n is unknown, three simultaneous equations can be formed.

Under the Casson model:

σ^(0.5)=σ₀ ^(0.5)+(PV·γ)^(0.5),  (29)

where σ is the shear stress, σ is the shear rate, σ₀ is the yield stress and PV is the plastic viscosity. Applying the Cox-Merz rule as before leads to the following simultaneous equations (assuming power law index n is known):

σ₀ ^(0.5)=(V _(L)(ω1)·ω1)^(0.5)−(PV·ω1)^(0.5),

σ₀ ^(0.5)=(V _(L)(ω2)·ω2)^(0.5)−(PV·ω2)^(0.5),  (30)

which can be solved to give estimates of the yield stress and plastic viscosity.

Other models of viscosity are known and this approach to finding the model parameters such as yield stress via viscosity measurements at multiple frequencies is not limited to the Herschel-Bulkley and Casson model examples presented here.

The yield stress can alternatively be obtained, using equations (18), (28) or (30) for example, without requiring a viscometer capable of vibrating (e.g. resonating) at multiple frequencies. Instead, two transducers may be located relatively close to each other in the yield stress material to eliminate or minimize the effect of any spatial variation in properties of the yield stress material such as any spatial variation of temperature of spatial variation of composition (e.g. concentration). The transducers are configured to vibrate at their own particular frequencies that are different from one another.

A measurement at a third frequency of operation, higher still than the first and second frequencies of operation, such as at 5 kHz or more, can provide other advantages. In particular, the higher third frequency generates a shear wave in the yield stress material with an even shorter penetration depth than the shear wave of the second (i.e. middle) frequency. This frequency shear wave will encounter much the same liquid region as the second (i.e. middle) frequency shear wave, but at a higher shear rate. In this case, the second and third frequency shear waves will largely penetrate only the highly liquid boundary layer space without any dominating yield stress effect because neither shear wave penetrates the solid region to any great extent. Therefore these frequency measurements can be used to calculate power law parameters for the material using, equation (22) for example.

Estimating the yield stress of a yield stress material can also be performed with a single measurement frequency of vibration, whereby a long propagating wave responds to yield stress variations without requiring a second wave at a second frequency of vibration. The flow can be modulated to induce change between solid and liquid states. These changes can then be correlated with the response of a single-frequency vibrator viscometer.

The flow modulation or modification can be performed by pump speed modulation, flow diversion, constriction modulation by a valve, stop start flow, stirrer/agitator, local vibrator device, transducer displacement in the fluid, vessel displacement, eccentric stirring.

Local flow changes can be induced by modulating the spatial position of: the viscometer body orientation to the flow, the viscometer relative to nearby surface, a local flow obstructer, and the fluid container/pipe/vessel.

FIG. 15 shows flow modulation by viscometer body orientation to the flow. Depending on the orientation, a greater or lesser portion of the material adjacent to the transducer will be sheltered from the flow by the transducer, affecting the amount viscosity readings. By taking viscosity measurements when at more than one orientation, the yield stress may be estimated using equation (18) above and setting ω1=ω2.

Other approaches may estimate the yield stress using equation (28) or equation (30).

For example, a change in flow velocity at a viscometer will likely result in change in viscosity measurement. It has been found that the yield stress of the fluid affects the degree to which the viscosity changes as a result of the change in flow velocity.

Thus many empirical models can be fit to calibrate data accounting for the change in viscosity as a result of a change in flow velocity (as measured by, for example, a flow meter).

For example, it has been found that the yield stress can be estimated based on the change in viscosity resulting from a change in flow velocity by evaluating the following expression:

σ₀=(ΔV·K ₀·(U ₁ ^(−K) ¹ −U ₂ ^(−K) ¹ )⁻¹)^(K) ² ,  (30A)

where U₁ and U₂ are the flow velocities, ΔV is the change in viscosity arising from the change in flow velocity, and K₀, K₁, and K₂ are empirical constants, which can be found by a standard calibration process.

The change in flow velocity is not limited to temporal changes in flow velocity. For example, the flow velocity can be a changed as a result of a change in location of the viscometer. Alternatively, two viscometers can be located in different locations. The two viscometers may encounter different flow velocities, for example, as a result of different pipe widths. The expression of equation (30A) may still be applicable and provide an estimate of the yield stress. Such an embodiment may be of particular advantage when multiple viscometers are available but limited to the same single resonant mode (resonant frequency). Such a program may have advantages in terms of economy of manufacture, particularly for MEMS and NEMS applications.

The invention is of course not limited to such a precise formula and other formula may be readily determined based on the principle of changing the flow around a viscometer between frequency and viscosity measurements.

FIG. 16 shows flow modulation by viscometer position relative to a nearby surface. Either or both of the surface and transducer can move. The distance from the boundary layer at the nearby surface affects the boundary layer around the transducer. By modulating the distance and continuing to take viscosity measurements, the yield stress may be determined.

FIG. 17 shows flow modulation by a local flow obstructer. Upstream of the transducer an obstruction changes position, which may include a change of orientation, relative to the transducer causing varying amounts of sheltering on the transducer. By modulating the position of the obstruction and continuing to take viscosity measurements, the yield stress may be determined. The movement of other obstructions may also affect the flow around the transducer. For example, movement of an obstruction downstream of the transducer may also affect flow around the transducer; what is important is that the obstruction has an effect on the flow around the transducer. For example, the obstruction may be in a vicinity of the transducer.

FIG. 18 shows flow modulation by modulating a container position, the container being the container in which the yield stress material is held along with the transducer.

Another approach to achieving flow modulation may be through vortex shedding, either by the transducer itself or by an obstruction upstream of the transducer. Vortex shedding provides an oscillating flow whereby vortices are created on the downstream side of the vortex shedding body, the vortices detaching periodically from either side of the body.

In these cases, the alteration of flow gives rise to a change of boundary layer formation and therefore presents different solid/liquid areas to the viscometer measurement region. This causes the measured viscosity to vary accordingly and the variation of viscosity can be scaled to yield stress. For example, in the case of a fluid with no yield stress, modulation of the flow will make little or no difference to the measured viscosity

Liquefaction may also have application in the determining of a yield stress or other fluid property of a yield stress fluid. Some yield stress fluids, and particularly those having the form of a granular suspension, can reduce their apparent yield stress when exposed to vibrations at a high enough frequency and a high enough amplitude. By this principle, air is commonly removed from wet concrete. The wet concrete is vibrated to cause a reduction in its yield stress allowing air bubbles trapped in the semi-solid matrix of the wet concrete to be released. Borrowing terminology from seismology, the inventors describe this process in yield stress materials as ‘liquefaction’.

An investigation into the liquefying effect of sufficiently strong oscillations in a yield stress material that may be helpful in understanding the present disclosure may be found in Balmforth et al “The Viscoplastic Stokes Layer” (Journal of Non-Newtonian Fluid Mechanics, volume 158, 2009, doi:10.1016/j.jnnfm.2008.07.008). Other material that may assist may be found in Vavreck “Flow Of Molten Milk Chocolate From An Efflux Viscometer Under Vibration At Various Frequencies And Displacements” (International Journal of Food Science & Technology, volume 39, 2004, doi:10.1111/j.1365-2621.2004.00805.x).

Vibration sources of an appropriate amplitude and frequency may be applied locally to the yield stress material surrounding the viscometer transducer to create liquefaction. This may be performed in flowing or static conditions.

Yield stress change is detected by either the reduction in the apparent viscosity or a change in frequency due to a change in mass loading at the transducer surface.

The vibration source may be the resonance of the transducer itself as used to make a viscosity measurement and in the context of the present disclosure, this is termed ‘intrinsic stimulation’. Alternatively it may be a separate vibration.

The combination of high frequency and high amplitude of vibration may cause local liquefaction at and around the transducer surface, detected as a viscosity and frequency change. Modulation of transducer vibration may be managed by control electronics, which may also allow the simultaneous measurement of viscous loss and resonator frequency. The method can be used with either a single transducer operating at one frequency, multiple transducers with different operating frequencies, or a multi-frequency resonator unit, i.e. a vibratory transducer having a plurality of resonant modes and being capable of resonating at multiple frequencies.

Alternatively, or additionally, the vibration to create the liquefaction can be a different mode than the vibration to make the viscosity measurement. For example, it can be a second torsional mode, or it can be a lateral or longitudinal vibrational mode of the transducer.

FIG. 19 illustrates liquefaction by intrinsic vibration. If the yield stress material is stationary, then there is no flow, and the only liquefied region of the yield stress material is that produced by the vibrations of the transducer. If the yield stress material is flowing around the transducer, then a boundary layer develops naturally as a result of the velocity gradients generated by the presence of the transducer. A portion of the liquefied region around the transducer will be due to the boundary layer and a portion of the liquefied region around the transducer will be due to the liquefaction.

Alternatively, or additionally, a liquefaction vibration source may be external to the vibrating viscometer (or viscometers in the case where multiple viscometers are used). In the context of the present disclosure, this is termed ‘extrinsic stimulation’.

The liquefaction vibration source may comprise one or more mechanical vibration devices connected to the vessel or one or more acoustic vibration sources mounted on or near the transducer or transducers, at a location offset from the transducer or transducers within the vessel, at the vessel wall or, where a free surface is present such as in the case of open-channel flow, above the free surface of the fluid.

FIG. 20 illustrates liquefaction by extrinsic vibration, whereby vibration to cause liquefaction is applied to a wall of the vessel or conduit in which the yield stress material is located.

FIG. 21 illustrates liquefaction by extrinsic vibration, whereby vibration to cause liquefaction is provided by vibration sources located elsewhere in the yield stress material.

As with intrinsic stimulation, the combination of a relatively high frequency and amplitude of vibration by the liquefaction vibration source causes liquefaction near the transducer surface and is detected as a viscosity and frequency change. Modulation of transducer vibration may be managed by control electronics according to the desired protocol of use, such as with a single transducer at a single frequency, multiple transducers, or a multi-frequency unit.

It is noted that extrinsic stimulation has previously been combined with conventional rotational viscometers for static materials in captive samples but not with a resonant viscometer. A particular advantage is provided by the use of a resonant viscometer with which the frequency change due to different mass loading of the resonant transducer has particular significance.

It is noted that the liquefaction techniques described above are particularly effective with yield stress materials that are granular suspensions rather than gel-like materials. Further background on such materials may be found in Hanotin et al “Viscoelasticity Of Vibrated Granular Suspensions” (Journal of Rheology, volume 59, 2015, doi:10.1122/1.4904421).

Moreover, while the above discussion makes particular reference to the viscosity being measured and used with the liquefaction techniques to obtain estimates of the properties of the yield stress material, it is noted that the techniques do not require an exact calculation or determination of the viscosity. Instead, the skilled person will recognize, based on the teaching in this disclosure, that other quantities can be used that vary with viscosity, such as the damping factor, or the Q-factor, or the degree of attenuation of vibrations between successive cycles or vibrational half-life. Others may be possible and may be within the scope of this disclosure. More generally, in this disclosure the term ‘degree of damping’ is used to represent the broad class of damping quantities including viscosity that have viscosity dependency and may be used in the techniques described herein.

Knowledge of the boundary layer may be used to provide improved techniques and apparatus for estimating fluid properties such as the yield stress, such as by manipulating the flow to change the boundary layer in the area around a viscosity transducer.

A boundary layer will usually form where a yield stress fluid with a velocity makes contact with a surface. The presentation of any surface to a flowing yield stress fluid dictates the formation of solid and liquid regions in the fluid. Promoting regions of substantial flow will result in higher velocity gradients and greater proportion of regions in which the shear stresses are greater than the yield stress, and thus result in liquefied regions. Where there is low or zero flow, the velocity gradients in the fluid will be low or zero, leading to regions where the shear stresses are below the yield stress and therefore solid regions.

Manipulation of flow may favourably create high shear regions or low shear regions—known as ‘dead zones’—in the path of propagating waves. This may improve the viscous and mass loading response of the resonant viscometer(s) to changes in fluid rheology. Solid regions close to the transducer will present the material better for liquefaction while also amplifying the influence of yield stress on the measurements of the low frequency viscometer.

Local flow modification can arise from recesses or protrusions such as ridges in the transducer or container geometry or by introducing obstructions or diversions to create regions sheltered from the flow—‘flow shadows’—leading to areas of low or zero shear stress and the formation of solid regions.

By creating contours or by altering the relative proximity of surfaces it is possible to modify the flow to deliberately develop areas of ‘solid’, ‘transitional’ and ‘liquid’ material in the environment in which waves are propagated and thereby influence the yield stress measurement. For example, this may improve the sensitivity of the yield stress measurement.

In particular, by promoting the formation of solid regions at or close to the surface of the transducer, shear waves propagating from the surface of the transducer may encounter solid regions and transitional regions closer to the surface of the transducer, increasing sensitivity to the solid and transitional regions, from which the yield stress is determined.

Since the propagation depth of the shear wave decreases with increasing frequency, this may allow the transducer to vibrate at higher frequencies while still propagating through solid regions. This may be particularly advantageous in that a greater choice in transducers is provided if the sensitivity to solid regions can be ‘tuned’ in this way. Alternatively, it may be desirable to avoid certain frequency bands due to, for example, plant noise, but still obtain a sensitivity to solid regions equivalent or similar to a transducer operating in a restricted band.

Detecting the onset of flow is important in many applications. In many measurement environments it is desirous to be able to detect the onset of flow, i.e. the point at which static material begins to flow. One particularly important application may be the detection of the flow developing as a result of a slow leak in a system.

However, it can be difficult to achieve this using conventional flow measurement techniques which are typically not sensitive to flow at rates close to zero. This is particularly the case with yield stress materials, which may appear solid and unmoving but may be undergoing a very slow or creeping flow.

Without wishing to be bound by theory, it is believed that, since a yield stress fluid requires a liquid boundary layer in order to flow, an incipient liquid layer will form at the very outset and this liquid is manifestly detectable as a sudden change in viscosity from the vibratory viscosity transducer.

In particular, if the viscosity in a static sample of a yield stress material is monitored using a vibratory viscosity transducer, the measured viscosity will be observed to increase substantially from the value at static conditions as the yield stress material begins to flow, before decreasing to below the value at static conditions as flow rate increases beyond a barely perceptible creeping flow.

The technique is found to be sensitive, responding to the mere application of pressure to a yield stress fluid.

Following this principle, it is possible to determine when a yield stress material fluid is no longer static, in that it has begun to creep slowly, by detection of an increase in apparent viscosity (i.e. damping), using a vibratory viscometer in a yield stress material.

Thus there is provided a technique for determining movement that would otherwise be considered imperceptible.

FIG. 22 illustrates schematically an example apparatus for analysing a fluid using one or more techniques of this disclosure. The apparatus comprises a first resonant viscosity transducer 100 and a second resonant viscosity transducer 120 located in relatively close proximity to the first resonant viscosity transducer in a fluid sample 5. The fluid sample 5 in this instance is a fixed volume of generally stationary fluid in a chamber 6, the fluid having a free surface 7, and portions of the viscosity transducers piercing the free surface 7 from above.

While the chamber 6 is drawn in FIG. 22 having a closed top wall, as may be necessary if its contents are pressurized, the chamber may equally be open from above. The fluid sample may therefore be at atmospheric pressure. In such a configuration, the viscosity transducers may be located above the chamber 7 such that at least portions of the transducers extend into the upper opening of the chamber 7 and contact the fluid sample 5 instead.

The resonant viscosity transducers are of a type described in U.S. Pat. No. 6,450,013, in which the transducers include a vibrating element configured to oscillate in a torsional mode. The vibrating element is immersed in the fluid and the viscosity is determined by correlation with the damping experienced by the element, i.e. the Q factor. In particular, each transducer comprises a transducer mounting 10, a semi-rigid connection member 12, a shaft 14 and a sense element 16. The shaft 14 and the sense element 16 are driven to vibrate torsionally with an angular frequency ω. The sense element 16 and the shaft 14 and sense element 16 are formed, at least substantially and possibly entirely, of a metal material such as a stainless steel. The sense element 16 and shaft 14 both have a circular cross-section, i.e. they are circularly symmetrical about the axis of oscillatory rotation. An example transducer that may be suitable for determining the viscosity via vibration at a frequency is the XL7 model viscometer manufactured by Hydramotion Ltd of Malton, UK.

The contents of the chamber 6 are pressurized to 10 Bar relative to atmospheric pressure. The sensing element is exposed to the viscous effect of the fluid in the sample 5. Increasing viscosity of the fluid causes an increased damping of the vibration in the transducer, resulting in a measurable reduced vibrational efficiency of the system.

In this apparatus, the first viscosity transducer 100 is specially designed, via choice of stiffness and mass or moment of inertia of the resonant system, to have a low resonant frequency at 400 Hz, i.e. an angular frequency of approximately 2513 rad/s. The second viscosity transducer 120 is specially designed to have a higher resonant frequency at 1500 Hz, i.e. an angular frequency of approximately 9425 rad/s.

To determine the viscosity at the resonant frequencies of the first and second viscosity transducers, the ‘Q factor’ of the vibration can be determined. The Q factor is a dimensionless parameter that indicates the level of damping of a resonator, wherein the level of damping is a function of the viscosity. In particular, it indicates the degree to which a resonator is underdamped. On a plot of frequency response, a high Q factor provides a high and narrow peak at the resonant frequency whereas a low Q factor provides a low and wide peak. Due to the change in width of the peak with damping, the Q factor can be defined as the ratio of the resonant frequency to the resonant bandwidth:

$Q = \frac{\omega_{R}}{\Delta\omega}$

wherein ω_(R) is the resonant frequency in radians per second and Δω is the Full Width at Half Maximum (FWHM), the bandwidth over which the power of the vibration is greater than half of the maximum (or equivalently the amplitude of vibration is greater than the maximum amplitude at resonance divided by √2), i.e. the bandwidth between the 3 dB points. The fluid viscosity is inversely proportional to the square of the Q factor and any constant of proportionality needed to compute the value of the viscosity measurement can be obtained by calibration with reference fluids of known viscosity.

It should be noted that the measurement of viscosity at or corresponding to a frequency of vibration may comprise making amplitude measurements at more than one frequency to estimate the Q factor but a single viscosity measurement is obtained at a frequency corresponding to the or a resonant frequency. For example, the bandwidth can be determined based on the frequencies required to cause the amplitude to drop to a factor of 1/√2 of the maximum amplitude at resonance. As a non-limiting example, the frequencies required to cause the amplitude to drop to a factor of 1/√2 of the maximum amplitude at resonance may be determined by performing a frequency sweep around the resonant frequency, but the skilled reader will recognize that the 3 dB point frequencies can be identified by various other techniques.

Another approach to determining the Q factor is to measure the amplitude of vibration at a series of frequencies around the resonant frequency and fit a parabola by the method of least squares to the frequency and amplitude values (or logarithms thereof). The 3 dB points can then be obtained as solutions to a quadratic equation based on the parabola of best fit to the measurements.

Another approach to determining the Q factor is by logarithmic decrement. By ceasing to drive the transducer and measuring the decay of vibrations, the Q factor may be determined by monitoring time series of the vibrations and determining the natural logarithm of the ratio of two successive peaks, A₁ and A₂, by the following expression:

$Q = {\frac{1}{2}{\sqrt{1 + \left( \frac{2\pi}{\ln \left( \frac{A_{1}}{A_{2}} \right)} \right)^{2}}.}}$

The first and second viscosity transducers 100, 120 both provide viscosity measurements corresponding to their resonant frequencies. The first and second viscosity measurements V1, V2, and first and second angular frequencies ω₁, ω₂ are provided to a processing module (not shown) which processes these measurements to provide an estimate of the yield stress using equation (18) above, i.e. evaluating

$\sigma_{0} = {\alpha^{\prime}\left( {{V\; 1} - {V\; 2\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}}} \right)}$

using a known value of n, the power law index, and a previously obtained empirical scaling constant α′.

FIG. 23 illustrates schematically a further example apparatus for analysing a fluid using one or more techniques of this disclosure. The apparatus comprises a single resonant viscosity transducer 300 that has been specially designed to be able to operate at two different frequencies, 400 Hz and 1500 Hz, i.e. angular frequencies of approximately 2513 rad/s and 9425 rad/s, as with the apparatus shown in FIG. 22. Such a transducer may be obtained by choice of stiffness and mass or moment of inertia for each of the elements of the transducer such that it has at least two resonant modes corresponding to the desired frequencies. The fluid sample 5 in this instance is a fixed volume of fluid in a chamber 6, the fluid having a free surface 7, and a portion of the viscosity transducer 300 piercing the free surface 7 from above. In contrast with the example shown in FIG. 22, the chamber 6 is provided with a paddle stirrer 8 which continuously agitates the fluid sample 5 such that the fluid is flowing past the immersed portion of the viscosity transducer 300.

The first and second viscosity measurements V1, V2, and first and second frequencies ω₁, ω₂ are provided to a processing module 18 which processes these measurements by evaluating equation (18) above using a known value of the power law index to provide an estimate of the yield stress σ₀.

FIG. 24 illustrates schematically a further example apparatus for analysing a fluid using one or more techniques of this disclosure, in which multiple resonant viscosity transducers 400, 420, 440 each vibrate at increasingly higher frequencies corresponding to their respective resonant frequencies. It will be appreciated that the ordering of the transducers in increasing resonant frequency is not essential and they may be in any arbitrary order. In this instance, the resonant frequencies are 400 Hz, 1500 Hz and 5000 Hz, i.e. angular frequencies of approximately 2513 rad/s, 9425 rad/s and 31416 rad/s. The fluid sample 5 in this instance is a flowing in a conduit 7 at an upstream average speed of 1 m/s. The viscosity transducers 400, 420, 440 extend from above through the wall of the conduit 7, the portions extending through the wall of the conduit 7 being in contact with the fluid sample as it flows past the viscosity transducers 400, 420, 440.

The first, second and third viscosity measurements V1, V2, V3 and first and second frequencies ω₁, ω₂, ω₃ are provided to a processing module (not shown) which processes these measurements to provide estimates of one or more fluid properties.

The processing module evaluates the power law index n from the measured V2, V3, ω₂, and ω₃ by substitution into equation (22) above, i.e.:

$n = {\frac{\log \left( {V\; 2\text{/}V\; 3} \right)}{\log \left( {{\omega 2}\text{/}{\omega 3}} \right)} + 1}$

Once the power law index is known, the yield stress can be evaluated using equation (18) above and V1, V2, ω₁, and ω₂, i.e.

$\sigma_{0} = {\alpha^{\prime}\left( {{V\; 1} - {V\; 2\left( \frac{\omega 1}{\omega 2} \right)^{n - 1}}} \right)}$

FIG. 25 illustrates schematically a further example apparatus for analysing a fluid using one or more techniques of this disclosure. The apparatus includes a single resonant viscosity transducer 500 that has been specifically designed to be able to operate at multiple frequencies: 400 Hz, 1500 Hz and 5000 Hz as per the embodiment shown in FIG. 24, i.e. angular frequencies of approximately 2513 rad/s, 9425 rad/s and 31416 rad/s. Such a transducer may be obtained by choice of stiffness and mass or moment of inertia for each of the elements of the transducer such that it has at least two resonant modes corresponding to the desired frequencies. The viscosity transducer 500 is mounted horizontally through the side wall 9 of a chamber and the portion extending through the side wall is immersed in the fluid sample 5. The viscosity transducer 500 has a particularly long shaft 14 such that the sense element 16 is spaced a particularly long distance from the walls 9 of the chamber.

The first, second and third viscosity measurements V1, V2, V3, and first and second frequencies ω₁, ω₂, ω₃ are provided to a processing module (not shown) which processes these measurements to provide estimates of one or more fluid properties.

As with the embodiment shown in FIG. 24, the processing module evaluates the power law index n and yield stress σ₀ from the measured V1, V2, V3, ω₁, ω₂, and ω₃ by substituting these values into equations (22) and (18) in turn.

FIG. 26 illustrates a variant apparatus for analysing a fluid using one or more techniques of this disclosure, in which multiple resonant viscosity transducers 400, 420, 440 each vibrate at increasingly higher frequencies corresponding to their respective resonant frequencies. The apparatus includes a single resonant viscosity transducer 600 that has been specifically designed to be able to operate at two frequencies: 400 Hz and 1500 Hz as per the second embodiment shown in FIG. 23, i.e. angular frequencies of approximately 2513 rad/s and 9425 rad/s. Such a transducer may be obtained by choice of stiffness and mass or moment of inertia for each of the elements of the transducer such that it has at least two resonant modes corresponding to the desired frequencies.

The viscosity transducer 600 comprises a sense element 16 on a shaft 14 extending from a threaded transducer mount 10. A processing module 18 is shown on the other side of the threaded transducer mount 10. The threaded portion of the transducer mount 10 engages with corresponding threads of a conduit 8. The sense element 16 and shaft extend along the axis of the conduit 8 facing upstream. The fluid-contacting portions of the viscosity transducer are formed of stainless steel (type 316) to resist corrosion and avoid contamination of the fluid. The conduit 8 includes a corner such that the fluid is diverted by the conduit walls to one side after it has passed the sense element 16 of the transducer 600, i.e. the viscosity transducer 600 is installed in an elbow section of the conduit 8.

The processing module 18 outputs the yield stress σ₀, and the fluid temperature. The fluid temperature is obtained from a temperature sensor included in the sense element 16 of the viscosity transducer 600. The yield stress σ₀ is obtained from the measured V1, V2, ω₁, and ω₂ by substitution into equation (18) as above.

In the examples shown in FIGS. 22 to 26, the fluid is variously at rest, in motion in a chamber, stirred by a paddle, and flowing in a conduit. It will be recognized that all of the example apparatus may be used in any of such fluid environments. For example, while FIG. 24 illustrates fluid flowing in a conduit past three separate viscosity transducers, these could be replaced by the single multimode viscosity transducer in FIG. 25 or any of the viscosity transducer arrangements in FIGS. 22 and 23, or some combination of multimode and single mode transducers.

The variant apparatus of FIG. 26 could equally have been included in any of the configurations of FIGS. 22 to 25, where the pipe fitting may be replaced with a more suitable means of attachment for extending through a bulkhead or conduit wall.

Alternatively, the above-described techniques may be implemented using another type of viscosity transducer, such as a vibrating-tube-type viscosity transducer. An example of a suitable transducer of this type is described in WO 2017001861 A1.

FIG. 27A shows a viscosity transducer 700 of the shaft-and-bob type comprising a shaft 714 and a bob 716 extending from a conduit wall 718 into the path of flowing yield stress fluid 5. The viscosity transducer 700 includes a pivot means 712 and is capable varying the orientation of the shaft 714 and bob 716 in the fluid 5.

In FIG. 27A, the shaft is aligned perpendicular to the direction of flow. In this configuration, a viscosity V1 and an angular frequency ω1 is obtained by torsional vibration of the shaft 714 and bob 716.

In FIG. 27B, the shaft 714 and bob 716 are rotated about pivot means 712 into the flow by an angle θ and a viscosity measurement is made at this orientation, leading to a viscosity V2 and an angular frequency ω2.

From these measurements, the yield stress of the fluid may be estimated using one of the approaches described in this disclosure, such as one of equations (18), (28) and (30).

Liquefaction may be introduced to liquefy regions of the yield stress material through the application of vibrations of sufficient amplitude and frequency.

In an embodiment, a single-frequency vibratory transducer is driven at its resonant frequency and the viscosity signal derived from shear losses is measured. Liquefaction is frequency-dependent and favours higher frequency so the frequency of vibration in this case may be of the order of 2 kHz but the invention is not limited to this frequency.

The liquefying vibration can be the actual vibration used for measurement, or any other torsional, lateral or longitudinal mode of vibration of the viscometer or the viscometer body.

In the case where the vibration is also the measurement mode, the amplitude A of vibration is modulated from a very low level, sufficient for measurement but not liquefaction and the value then recorded. Where there material has yield stress, increasing the amplitude will lead to local liquefaction nearest the surface of the resonator which may be detected in two ways.

In a first approach, a change in measured viscosity as the liquid layer is formed, ΔV_(AMP). Typically this is a reduction in viscosity as the apparent high viscosity of a semi-solid matrix is broken down. In other cases, where the yield stress material has a more gel-like matrix, the impedance to the measuring wave may increase as a lossy liquid is formed and the viscous loss may actually increase. The viscous loss change is then scaled (k_(LV)) against wave (resonator) amplitude to provide an estimate of yield stress.

A viscometer is driven at a first amplitude A_(A), sufficient to liquefy the yield stress material, leading to measurements of viscosity V_(A) and resonant frequency F_(A). The viscometer is then driven at a second, lower, amplitude A_(B), leading to measurements of viscosity V_(B) and resonant frequency F_(B). This leads to ‘delta’ values of ΔA=A_(A)-A_(B), ΔV_(AMP)=V_(A)−V_(B), and ΔF=F_(A)−F_(B), representing the step change in amplitude and the resulting changes in viscosity and frequency. An empirical estimate of the yield stress can be made by way of the following expression:

σ₀ =ΔV _(AMP) k _(LV) /ΔA  (31)

In a second approach, that may be implemented as an alternative or in addition to the first (i.e. viscosity change) approach, there may be an increase in frequency F as the more liquefied layer presents a lower mass loading m_(F) than the solid structure

$\begin{matrix} {F = {\frac{1}{2\pi}\sqrt{K_{v}\text{/}\left( {m_{0} + m_{F}} \right)}}} & (32) \end{matrix}$

in which K_(v) and m₀ are stiffness and mass parameters of the mechanical vibrating system.

The frequency change ΔF_(AMP) is then scaled (k_(LF)) against the wave (resonator) amplitude to provide an estimate of yield stress by the following expression.

σ_(p) =ΔF _(AMP) k _(LF) /ΔA  (33)

In an embodiment, vibration at a second (lower) frequency is additionally employed which does not contribute to liquefaction. This provides a reference signal to ratiometrically reduce or eliminate systematic errors caused by temperature or changes in fluid viscosity. By dividing by the viscosity or frequency measurements corresponding to a lower frequency (which has negligible liquefaction effect) errors which might be caused by the natural common mode changes of temperature or viscosity in the fluid are reduced or eliminated.

In particular, a vibratory viscometer is driven at a first resonant mode at a first amplitude A_(A), sufficient to liquefy the yield stress material, leading to measurements of viscosity V_(A-HIGH) and resonant frequency F_(A-HIGH). The viscometer is then driven at the first resonant mode at a second, lower, amplitude A_(B), leading to measurements of viscosity V_(B-HIGH) and resonant frequency F_(B-HIGH). This leads to a ‘delta’ value of ΔA=A_(A)−A_(B), representing the step change in amplitude.

A second vibratory viscometer is also driven at a at a second resonant mode (corresponding to a resonant frequency lower than at the first resonant mode) at an amplitude insufficient to cause liquefaction (e.g. A_(A)), leading to measurements of viscosity V_(A-LOW) and/or resonant frequency F_(A-LOW). The viscometer is then driven at the first resonant mode at an amplitude insufficient to cause liquefaction (e.g. A_(A)), leading to measurements of viscosity V_(B-LOW) and/or resonant frequency F_(B-LOW).

Alternatively, a single multi-mode viscometer can be used to take all of the measurements.

In the case of viscosity, the viscosity ratio change ΔV_(RL), may then be scaled (k_(LRL)) against change in drive amplitude to provide an improved estimate of yield stress in which common mode errors, such as errors due to temperature change, are reduced or eliminated:

ΔV _(RL) =V _(A-HIGH) /V _(A-LOW) −V _(B-HIGH) /V _(B-LOW)

σ₀ =ΔV _(RL) ·k _(LRL) /ΔA  (34)

Alternatively or additionally, for frequency, the frequency ratio change ΔF_(RL) may then be scaled (k_(LRF)) against change in drive amplitude to provide an improved estimate of yield stress:

ΔF _(RL) =F _(A-HIGH) /F _(A-LOW) −F _(B-HIGH) /F _(B-LOW)

σ₀ =ΔF _(RL) k _(LRF) /ΔA  (35)

It is noted that, according to the techniques of the present disclosure, the frequency may be measured as an intermediate or accompanying step in the process of measuring the viscosity of the yield stress material. For example, identifying a Q-factor may require the resonant frequency to be identified. On the other hand, the resonant frequencies might be determined without also measuring viscosities.

Alternatively, it is not required that the viscometers (transducers) vibrate at resonance, particularly in the case where a change in viscosity is used to estimate the yield stress. For example, a viscometer may operate at a frequency or frequencies away from a resonant mode, such as a lower frequency, and obtain measurements of viscosity. In such cases away from resonance, frequency-based Q-factor approaches such as using bandwidth to determine a Q-factor and from the Q-factor a viscosity, may not be appropriate. It may be more appropriate to determine viscosity based on time-series approaches such as the log-decrement method or by a consideration of drag forces.

In such cases, the yield stress may be estimated by vibrating a vibratory transducer in the yield stress fluid at a first frequency and making a first measurement of the viscosity; providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrating a vibratory transducer the first frequency and making a second measurement of the viscosity; estimating the yield stress of the yield stress fluid based on the first and second measurements of viscosity. The vibration to liquefy yield stress fluid around the one or more transducers is provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement. The yield stress may be estimated based on the difference between the first and second measurements of the viscosity scaled by the difference in amplitude between the first and second measurements such as by equation (31).

In another embodiment, a vibratory transducer is vibrated at a second frequency in the yield stress fluid and a third measurement is made of the viscosity, the second frequency being lower than the first frequency, and a vibratory transducer is vibrated at the second frequency in the yield stress fluid and a fourth measurement is made of the viscosity, wherein the yield stress is estimated based on difference between the first viscosity measurement scaled by the third viscosity measurement and the second viscosity measurement scaled by the fourth viscosity measurement, scaled by the difference in amplitude between the first and second measurements.

It may not be necessary to take separate third and fourth measurements as they may be the same—it is believed that amplitude of vibration has little effect on viscosity or resonant frequency provided the amplitude (or amplitude-frequency product) is sufficiently low that liquefaction does not occur.

For improved reduction in temporal fluctuations, the third and fourth measurements of viscosity or frequency may be made simultaneously with the first and second measurements respectively, in the case of a multi-frequency (e.g. multi-mode) viscometer, or immediately before or immediately after.

The invention is not limited to step changes in amplitude. In other embodiments, one or more step variations may be employed, including a series of gradual variations such as with an amplitude sweep, or periodic variation of amplitude according to sine waves, triangle waves, sawtooth waves, square waves or other periodic (or even aperiodic) signals. These changes may be detected as viscosity and/or frequency changes accordingly.

Moreover, as noted previously, these techniques are not limited to determining the viscosity per se. Other viscosity-like quantities may be used, such as damping factor, Q-factor, attenuation rate, that are associated with the viscosity and vary with it due to strong viscosity dependence. These quantities may be referred to more generally, within the context of this disclosure, as representing the degree of damping.

The source of vibration to liquefy the yield stress fluid around the viscometer is not required to originate from the viscometer.

FIG. 20 shows external mechanical vibration of structure coupled to fluid in any of torsional, lateral and longitudinal modes. Independent to the viscometer, the vibration may be provided by mechanical stimulation of the container or infrastructure surrounding the fluid as shown in FIG. 20.

Alternatively, or additionally, fluid-borne acoustic vibration can be the liquefying source as shown in FIG. 21.

In all cases the amplitude and frequency of these source can be controlled and modulated and the same measurements made as described above in respect of yield stress measurement from intrinsic liquefaction.

Alternatively, the above-described techniques are not limited to a particular type of vibrating transducer and may be implemented using a shaft-and-bob-type transducer or another type of viscosity transducer, such as a vibrating-tube-type viscosity transducer. An example of a suitable transducer of this type is described in WO 2017001861 A1.

FIGS. 28 to 33 show variations of this approach with viscosity transducers of the bob-and-shaft type, with the shaft axially aligned with the direction of flow, the bob facing upstream.

FIG. 28 shows a transducer with a smooth profile in open flow of yield stress material. The profile of the bob is rectangular in cross-section (i.e. cylindrical due to symmetry about axis), its width (i.e. diameter) reducing steadily to a point at its conical upstream end and its width (i.e. diameter) reducing steadily to the narrower width of the shaft at its frustro-conical downstream end. A boundary layer develops around the transducer, and thus waves propagating in the yield stress material from the transducer enter a liquid region followed by a solid region.

FIG. 29 shows the transducer of FIG. 28 in a flow of yield stress material through a pipe or vessel, the walls of the pipe or vessel being axially aligned with the flow direction and the shaft of the transducer, wherein a further boundary layer develops at the walls of the pipe or vessel. Therefore waves propagating in the yield stress material from the transducer enter a liquid region followed by a solid region followed by a further liquid region at the walls of the pipe or vessel.

FIG. 30 shows a transducer with a contoured profile in an open flow of yield stress material. The transducer is different from the transducer of FIG. 28 in that the rectangular (i.e. cylindrical) portion of the bob is provided with recesses, which are circumferential grooves in this case. Material in the recesses is sheltered from the flow and low shear stresses are generated within the recesses, with the result that the material in the recesses remains solid. A liquid boundary layer develops around the transducer but does not extend into the recesses. Therefore waves propagating in the yield stress material from the transducer enter a solid region in the recesses followed by a liquid region in the boundary layer at the transducer away from the recesses, followed by a further solid region in the open flowing yield stress material.

FIG. 31 shows the transducer of FIG. 30 with a contoured profile in a flow of yield stress material through a pipe or vessel, wherein a further boundary layer develops at the walls of the pipe or vessel. Therefore waves propagating in the yield stress material from the transducer enter a solid region in the recesses of the transducer followed by a liquid region in the boundary layer at the transducer away from the recesses followed by a solid region followed by a further liquid region at the walls of the pipe or vessel.

FIG. 32 shows the smooth transducer of FIG. 28 in a flow of yield stress material through a pipe or vessel, the walls of the pipe or vessel being axially aligned with the flow direction and the shaft of the transducer, wherein recesses are provided in the walls of the pipe or vessel, in this case in the form of circumferential grooves. Material in the recesses is sheltered from the flow and low shear stresses are generated within the recesses, with the result that the material in the recesses remains solid. A liquid boundary layer develops at the walls of the pipe or vessel at but does not extend into the recesses. A liquid boundary layer also develops around the transducer. Therefore waves propagating in the yield stress material from the transducer enter a liquid region followed by a solid region followed by another liquid region at the walls followed by another solid region in the recesses at the walls.

FIG. 33 shows the transducer of FIG. 30 with a contoured profile in a flow of yield stress material through a pipe or vessel, wherein a further boundary layer develops at the walls of the pipe or vessel, the walls of the pipe or vessel being axially aligned with the flow direction and the shaft of the transducer, wherein recesses are provided in the walls of the pipe or vessel, in this case in the form of circumferential grooves. Therefore waves propagating in the yield stress material from the transducer enter a solid region in the recesses of the transducer followed by a liquid region in the boundary layer at the transducer away from the recesses followed by a solid region followed by a further liquid region at the walls of the pipe or vessel followed by another solid region in the recesses at the walls.

FIG. 34 shows the transducer of FIG. 28 in open flow of yield stress material but, unlike in FIG. 28, the transducer is not axially aligned with the flow direction. It is instead aligned perpendicular to the flow direction. In this case, the liquid boundary layer develops on the upstream side of the transducer but the downstream side of the transducer is sheltered from the flow and may remain solid. Thus different propagation paths will be experienced by waves propagating from different sides of the transducer. The orienting of a transducer perpendicular to the flow direction in this way may be advantageously combined with any other techniques in this disclosure, including, but not limited to, the estimating of yield stress via difference in viscosities (e.g. by equation (18)), liquefaction).

FIG. 35 shows the transducer of FIG. 34 aligned perpendicular to flow in a pipe, wherein the transducer is partially retracted into a recess of the pipe. In this case, the liquid boundary layer develops on the upstream side of the transducer but the downstream side of the transducer is sheltered from the flow and may remain solid. In addition, a liquid boundary layer develops at the pipe wall upstream of the transducer. Again, different propagation paths will be experienced by waves propagating from different sides of the transducer, with the upstream propagating waves potentially experiencing the effect of the liquid boundary layer at the transducer and at the wall. In this case, the shaft is aligned perpendicular to the flow direction but other embodiments may have the shaft aligned with (parallel to) the flow direction or at some angle intermediate between perpendicular and parallel.

FIG. 36 shows a transducer aligned perpendicular to flow in a pipe wherein the transducer is a different shaft-and-bob-type transducer than FIG. 35. The bob is in the form of a disc or short cylinder aligned axially with the shaft. As with FIG. 35, the transducer is partially retracted into a recess of the pipe wall and the disc-like bob is located within what would ordinarily have been the boundary layer at the pipe wall.

FIG. 37A-E shows a series of five bob profiles as possible designs for a transducer of the shaft-and-bob type.

FIG. 37A shows a bob profile that is smooth sided and similar to the transducer of FIG. 22.

FIG. 37B shows a bob profile that has rectangular-shaped recesses on either side, representing circumferential grooves of rectangular profile.

FIG. 37C shows a bob profile that has quadrilateral-shaped recesses on either side that are widest at the outermost edge and narrow toward the axis, representing a circumferential groove with a flat upstream edge aligned transverse to the axis making right-angled corner with the bottom of the recess on the upstream side of the groove and the groove widening from the bottom of the recess in a frustro-conical fashion to the outer diameter of the bob at the downstream side.

FIG. 37D shows a bob profile that has triangle-shaped recesses on either side that are widest at the outermost edge and narrow toward the axis, representing a circumferential groove with a flat upstream edge aligned transverse to the axis, widening from the bottom of the recess in a frustro-conical fashion to the outer diameter of the bob at the downstream side.

FIG. 37E shows a bob profile that has recesses where the local diameter smoothly varies from a maximum to a minimum in a serpentine manner similar to a sinusoidal edge profile.

FIG. 38 shows two further bob designs for a transducer, in which the recesses/grooves take the form of an abrupt change in diameter. In the uppermost bob profile of FIG. 38, the diameter decreases sharply at a location along the axis as the flow moves downstream. This creates a recessed region in which flowing material is sheltered downstream of the diameter change and so the material may remain solid in this region. In the lowermost bob profile of FIG. 38, the diameter increases sharply at a location along the axis as the flow moves downstream. As the material flows around the transducer, material immediately upstream of the diameter change is retained in a recess and remains solid. Thus in FIG. 38, the recess may be described as relative to a streamline or flow path around the bob, or the bob profile may be viewed as having a ridge, before or after which yield stress material may be retained.

Alternatively, or additionally, the bob may be provided with ridges or recess extending to some axial extent along the bob. This may provide advantages because, for rotational vibrations, a slip layer or surface may form at the resonator surface, wherein the resonator is lubricated by the separation of phases in the yield stress material. If the shear region can be transferred further out from the resonator, then problems due to slip at the resonator surface are reduced or eliminated.

The retention of solid matter at the surface reduces or eliminates the formation of a slip layer at the sensor surface and thus provides improved coupling of the sensor to the fluid.

FIGS. 39A and 39B show bob shapes in a perspective view wherein smooth-bottomed grooves of semi-circular profile are provided spaced around the circumference of the bob extending axially along the bob.

FIG. 39A shows a variant, termed in this disclosure a ‘rifled’ bob, in which the path of the grooves has an axial and tangential component, leading to a helical path around the outer surface of the bob.

FIG. 39B shows another variant in which the grooves of semi-circular profile are axially aligned with the bob.

FIG. 39C shows a cross-section transverse to the axis of the bob that is applicable to either of the bobs shown in FIGS. 33A and 33B, showing solid yield stress material filling the grooves. The slip region is being moved to a circle taking in the outer extent of the bob, wherein the slip region is predominantly in in the yield stress material.

The ‘rifled’ bob has a further advantage in that, when axially aligned in the flow direction, it also may shelter material from the flow in its helical grooves in a similar manner to the bobs shown in FIG. 36B-E.

FIGS. 40A and 40B show a further bob variant in cross-sectional and perspective views respectively, in which the shear region is moved away from the resonator surface, in this case by providing vanes in the form of rectangular plates extending outward from the bob, the axis of the bob being coplanar with the vanes. Depending on the size of the vanes, the shear region can be moved as far as desired from the resonator surface, providing more flexibility in design than with grooves, which are limited by the dimensions of the bob.

Alternatively, the above-described techniques may be implemented using another type of viscosity transducer, such as a vibrating-tube-type viscosity transducer. An example of a suitable transducer of this type is described in WO 2017001861 A1. The tube internal surface may be provided with vanes, ridges, or recesses. aligned axially, tangentially, or helically around the tube.

In a further embodiment, a viscosity transducer (such as the viscosity transducer 100 in FIG. 22) is used to detect the beginning of flow in an initially static yield stress material.

A series of viscosity (i.e. damping) readings are taken using the viscosity transducer every 5 seconds for an ongoing sampling period.

Each successive viscosity measurement V_(n) is compared with its previous viscosity reading V_(n-1). If V_(n) is greater than V_(n-1) by more than a threshold, e.g. if V_(n)>ε·V_(n-1), where ε is a threshold ratio that is greater than 1 (e.g. 2), then the fluid is determined to have begun to flow.

In a further embodiment, two or more successive samples are averaged together, and the averaged samples are compared with preceding averages of successive samples, with the averaging effect serving as a low-pass filter to reduce the effect of noise and non-physical spikes and transients in the detection of flow.

Alternatively, the above-described techniques may be implemented using another type of viscosity transducer, such as a vibrating-tube-type viscosity transducer. An example of a suitable transducer of this type is described in WO 2017001861 A1.

In interpreting the disclosure, all terms should be interpreted in the broadest possible manner consistent with the context. In particular, the terms “comprises” and “comprising” should be interpreted as referring to elements, components, or steps in a non-exclusive manner, indicating that the referenced elements, components, or steps may be present, or utilized, or combined with other elements, components, or steps that are not expressly referenced. In the context of this disclosure, the term “based on” does not mean “based only on,” unless expressly specified otherwise. In other words, the term “based on” describes both “based only on” and “based at least on.” The term “determining” encompasses a wide variety of actions and, therefore, “determining” can include calculating, computing, processing, deriving, investigating, looking up (e.g., looking up in a table, a database or another data structure), ascertaining and the like. Also, “determining” can include receiving (e.g., receiving information), accessing (e.g., accessing data in a memory) and the like. Also, “determining” can include resolving, selecting, choosing, establishing and the like.

The methods, process and algorithms that have been described may be stored as one or more instructions on a processor-readable or computer-readable medium. The term “computer-readable medium” refers to any available medium that can be accessed by a computer or processor. By way of example, and not limitation, such a medium may comprise RAM, ROM, EEPROM, flash memory, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store desired program code in the form of instructions or data structures and that can be accessed by a computer. Disk and disc, as used herein, includes compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk and Blu-ray® disc where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. It should be noted that a computer-readable medium may be tangible and non-transitory. In the context of this disclosure, the term “code” may refer to software, instructions, code or data that is/are executable by a computing device or processor.

A processing module may comprise a computer including a processor for processing data and controlling systems according to the techniques and apparatus of the present disclosure. Alternatively, or additionally, the processing module may comprise electronic circuitry to perform the same functions, such as in the form of an FPGA (field programmable gate array) circuit and/or an ASIC (application specific integrated circuit), and/or a microcontroller.

Software or instructions or data may also be transmitted over a transmission medium. For example, if the software is transmitted from a website, server, or other remote source using a coaxial cable, fibre optic cable, twisted pair, digital subscriber line (DSL), or wireless technologies such as infrared, radio, and microwave, then the coaxial cable, fibre optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave are included in the definition of transmission medium.

The above detailed description of embodiments of the invention is not intended to be exhaustive or to limit the invention to the precise form disclosed above. While specific embodiments of, and examples for, the invention are described above for illustrative purposes, various equivalent modifications are possible within the scope of the invention, as those skilled in the relevant art will recognize. For example, while processes or blocks are presented in a given order, alternative embodiments may perform routines having steps, or employ systems having blocks, in a different order, and some processes or blocks may be deleted, moved, added, subdivided, combined, and/or modified. Each of these processes or blocks may be implemented in a variety of different ways. Also, while processes or blocks are at times shown as being performed in series, these processes or blocks may instead be performed in parallel, or may be performed at different times.

The teachings of the invention provided herein can be applied to other systems, not necessarily the system described above. The elements and acts of the various embodiments described above can be combined to provide further embodiments.

The headings provided herein are for convenience only and do not necessarily affect the scope or meaning of the aspects of this the disclosure defined by the claims.

Some embodiments have been described. These embodiments are presented by way of example only and are not intended to limit the scope of the disclosure. Indeed, the novel methods, apparatus and systems described herein may be embodied in a variety of other forms. It should be apparent to those skilled in the art that many more modifications besides those already described are possible without departing from the inventive concepts herein. For example, all methods described in the present disclosure may be alternatively embodied in apparatus for performing such methods, such as an apparatus comprising means for carrying out each step of such methods. As another example, all methods described in the present disclosure may be alternatively embodied in the form of a non-transitory (tangible) computer readable medium having instructions stored thereon that, when executed by a processor cause the processor to carry out the corresponding method.

While endeavouring in the foregoing specification to draw attention to those features of the invention believed to be of particular importance, it should be understood that the applicant claims protection in respect of any patentable feature or combination of features referred to herein, and/or shown in the drawings, whether or not particular emphasis has been placed thereon. 

1. A method of estimating the yield stress of a yield stress fluid using one or more vibratory transducers, the method comprising: vibrating a vibratory transducer in the yield stress fluid and making a first measurement; providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrating a vibratory transducer in the yield stress fluid and making a second measurement; estimating the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency.
 2. The method of claim 1, wherein the vibration to liquefy yield stress fluid around the one or more transducers is provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement.
 3. The method of claim 2, wherein the yield stress is estimated based on the difference between the first and second measurements, the difference scaled by the difference in amplitude between the first and second measurements.
 4. The method of claim 2, further comprising obtaining third and fourth measurements at a lower frequency and without liquefaction, wherein the yield stress is estimated based on the difference between the first measurement scaled by the third measurement and the second measurement scaled by the fourth measurement, said difference being scaled by the difference in amplitude between the first and second measurements.
 5. The method of claim 4, wherein the first, second, third and fourth measurements are made using a single multi-mode vibratory transducer.
 6. The method of claim 4, wherein the first and second measurements are made using a first vibratory transducer and the third and fourth measurements are made using a second vibratory transducer.
 7. The method of any of claims 1 to 6, wherein the yield stress fluid is static.
 8. The method of any of claims 1 to 6, wherein the yield stress fluid is flowing.
 9. The method of claim 1, wherein the yield stress fluid is in a vessel or conduit having a wall and the vibration to liquefy said portion of the yield stress fluid around the one or more vibratory transducers is provided by vibrating the wall of the vessel or conduit.
 10. The method of claim 1, wherein the yield stress fluid is in a vessel or conduit having a wall and the vibration to liquefy said portion of the yield stress fluid around the one or more vibratory transducers is provided by vibrating at least one or more additional transducers in the yield stress fluid.
 11. The method of any of claims 1 to 8, wherein the yield stress fluid is in a vessel or conduit having a wall and the vibration to liquefy said portion of the yield stress fluid around the one or more vibratory transducers is provided by vibrating at least one of the one or more vibratory transducers.
 12. The method of claim 11, wherein the vibration of at least one of the one or more vibratory transducers to liquefy said portion of the yield stress fluid around the one or more transducers is in the same vibratory mode as the vibratory mode used to obtain measurements.
 13. The method of claim 11, wherein the vibration of at least one of the one or more vibratory transducers to liquefy said portion of the yield stress fluid around the one or more transducers is in a different vibratory mode from the vibratory mode used to obtain measurements.
 14. The method of claim 13, wherein the vibration of at least one of the one or more vibratory transducers to liquefy said portion of the yield stress fluid around the one or more transducers is in a torsional vibratory mode and the vibratory mode used to obtain measurements is in a different torsional vibratory mode.
 15. The method of claim 13, wherein one of the modes of vibrations to liquefy said portion of the yield stress and to obtain measurements is a torsional mode and the other mode is a lateral or longitudinal mode.
 16. The method of any of claims 1 to 15, wherein the one or more vibratory transducers comprise a vibrating tube transducer.
 17. The method of any of claims 1 to 9 and 11 to 18, wherein the yield stress fluid is flowing relative to the one or more vibratory transducers and the one or more vibratory transducers comprise a vibratory transducer configured to vibrate in a torsional mode about an axis that is not aligned with the flow direction.
 18. The method of claim 17, wherein the axis is substantially perpendicular to the flow direction.
 19. The method of any of claims 1 to 18, wherein measurements of the degree of damping of the yield stress fluid comprise measurements of the viscosity of the yield stress fluid.
 20. An apparatus for estimating the yield stress of a yield stress fluid using one or more vibratory transducers, the apparatus comprising: one or more vibratory transducers; a processing module configured to: vibrate a vibratory transducer in the yield stress fluid and make a first measurement; provide a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrate a vibratory transducer in the yield stress fluid and make a second measurement; estimate the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency.
 21. The apparatus of claim 20, wherein the vibration to liquefy yield stress fluid around the one or more transducers is provided by the making of the second measurement at an increased amplitude of vibration relative to the first measurement.
 22. The apparatus of claim 21, wherein the yield stress is estimated based on the difference between the first and second measurements scaled by the difference in amplitude between the first and second measurements.
 23. The apparatus of claim 21, wherein the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the measurements are made of the resonant frequency, and the apparatus is further configured to make first and second measurements of the degree of damping based respectively on the first and second measurements of the resonant frequency, wherein the yield stress is estimated based on the difference between the first and second measurements of the degree of damping scaled by the difference in amplitude between the first and second measurements of resonant frequency.
 24. The apparatus of claim 21, wherein the processing module is further configured to obtain third and fourth measurements at a lower frequency and without liquefaction, wherein the yield stress is estimated based on the difference between the first measurement scaled by the third measurement and the second measurement scaled by the fourth measurement, said difference being scaled by the difference in amplitude between the first and second measurements.
 25. The apparatus of claim 24, comprising a multi-mode vibratory transducer and the first, second, third and fourth measurements are made using the multi-mode vibratory transducer.
 26. The apparatus of claim 24, comprising a first vibratory transducer and a second vibratory transducer, wherein the first and second measurements are made using the first vibratory transducer and the third and fourth measurements are made using the second vibratory transducer.
 27. The apparatus of any of claims 20 to 26, configured to estimate the yield stress of a static yield stress fluid.
 28. The apparatus of any of claims 20 to 26, configured to estimate the yield stress of a flowing yield stress fluid.
 29. The apparatus of any of claims 20, 27 and 28, wherein the yield stress fluid is in a vessel or conduit having a wall and the vibration to liquefy said portion of the yield stress fluid around the one or more vibratory transducers is provided by vibrating the wall of the vessel or conduit.
 30. The apparatus of any of claims 20, 27 and 28, wherein the yield stress fluid is in a vessel or conduit having a wall and the vibration to liquefy said portion of the yield stress fluid around the one or more vibratory transducers is provided by vibrating at least one or more additional transducers in the yield stress fluid.
 31. The apparatus of any of claims 20 to 28, wherein the yield stress fluid is in a vessel or conduit having a wall and the vibration to liquefy said portion of the yield stress fluid around the one or more vibratory transducers is provided by vibrating at least one of the one or more vibratory transducers.
 32. The apparatus of claim 31, wherein the vibration of at least one of the one or more vibratory transducers to liquefy said portion of the yield stress fluid around the one or more transducers is in the same vibratory mode as the vibratory mode used to obtain measurements.
 33. The apparatus of claim 31, wherein the vibration of at least one of the one or more vibratory transducers to liquefy said portion of the yield stress fluid around the one or more transducers is in a different vibratory mode from the vibratory mode used to obtain measurements.
 34. The apparatus of claim 33, wherein the vibration of at least one of the one or more vibratory transducers to liquefy said portion of the yield stress fluid around the one or more transducers is in a torsional vibratory mode and the vibratory mode used to obtain measurements is in a different torsional vibratory mode.
 35. The apparatus of claim 33, wherein one of the modes of vibrations to liquefy said portion of the yield stress and to obtain measurements is a torsional mode and the other mode is a lateral or longitudinal mode.
 36. The apparatus of any of claims 20 to 35, wherein the one or more vibratory transducers comprise a vibrating tube transducer.
 37. The apparatus of any of claims 20 to 35, wherein the yield stress fluid is flowing relative to the one or more vibratory transducers and the one or more vibratory transducers comprise a vibratory transducer configured to vibrate in a torsional mode.
 38. The apparatus of claim 37, wherein the torsional mode is about an axis that is not aligned with the flow direction.
 39. The apparatus of claim 38, wherein the axis is substantially perpendicular to the flow direction.
 40. An apparatus for estimating the yield stress of a yield stress fluid using one or more vibratory transducers, the method comprising: means for vibrating a vibratory transducer in the yield stress fluid and making a first measurement; means for providing a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; means for while said portion of the yield stress material is liquefied, vibrating a vibratory transducer in the yield stress fluid and making a second measurement; means for estimating the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency, wherein, optionally, the apparatus comprises means for performing a method according to any of claims 2 to
 19. 41. A non-transitory computer-readable medium having stored thereon instructions that, when executed by one or more processors, cause the one or more processors to: vibrate a vibratory transducer in the yield stress fluid and making a first measurement; provide a vibration to liquefy at least a portion of the yield stress fluid around the one or more vibratory transducers; while said portion of the yield stress material is liquefied, vibrate a vibratory transducer in the yield stress fluid and make a second measurement; estimate the yield stress of the yield stress fluid based on the first and second measurements, wherein, either: the vibratory transformer is vibrated in both cases at a first frequency and the measurements are made of the degree of damping in the yield stress fluid, or the vibratory transformer is vibrated in both cases at resonance in a first resonant mode and the first and second measurements are made of the resonant frequency, wherein, optionally, the instructions cause the one or more processors to perform a method according to any of claims 2 to
 19. 